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Claims: 



1 1 . A method for detecting anomalous phase measurements in a satellite differential 

2 navigation system in which a first satellite receiver and a second satellite receiver each compute 

3 phase measurements on a plurality of satellite channels, said method comprising the steps, 

4 performed for each of a plurality of iterations, of: 

5 a) generating a combined phase difference vector; 

6 b) generating a phase mismatch vector representing the difference between said 

7 combined phase difference vector and an estimated combined phase difference vector; 

8 c) generating an averaged estimate vector representing the averaged estimate of said 

9 phase mismatch vector over said channels; 

10 d) generating a residual vector representing the difference between said phase mismatch 

1 1 vector and said averaged estimate vector; 

12 e) generating a vector of controlling signals by linear transformation of said residual 

1 3 vector and said averaged estimate vector; 

14 f) generating said estimated combined phase difference vector by successively storing 

1 5 components of said vector of controlling signals for each of said channels; and 

16 g) detecting anomalous phase measurements by analyzing said residual vector. 



1 2. The method of claim 1 wherein said combined phase difference vector * c comprises 

2 phase differences (^ci..^Cj..^Cf^) for a j-th satellite channel, and wherein step a) of generating a 

3 combined phase difference vector further comprises the step of: 

4 calculating said phase difference vector f>cast>c = 4> -4> where vector 4> 

5 comprises the full phases measured by one of said satellite receivers for each j-th satellite 

6 channel (^bi -^ bj- -^bn) and vector 4> ^ comprises the full phases measured by the other satellite 

7 receiver for each j-th satellite channel i?/..* 

1 3. The method of claim 1 wherein said first and second satellite receivers are dual 

2 frequency (fl and f2) receivers and wherein said combined phase difference vector * c 

3 comprises phase differences (^*c/..*C;..*ca^) for a j-th satellite channel, and wherein step a) of 

4 generating a combined phase difference vector further comprises the step of: 

5 calculating said phase differences * cj of said phase difference vector * c as 
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6 *Cj= - - . *^j)(f2 / fl) 

7 where 4>®2j and represent the full phase measured by one of the satellite receivers at each j- 

8 th satellite channel at frequencies f2 and fl respectively, and 

9 where and represent the full phase measured by the other of the satellite receivers at 
10 each j-th satellite channel at frequencies f2 and fl respectively. 

1 4. The method of claim 1 wherein said step c) of generating an averaged estimate vector 

2 further comprises the step of: 

3 calculating said average estimate vector as 

4 N N 

5 A^ = ( 5:wj'A*Cj)/( Ewj) 
6 

7 where A*Cj represents the components of said phase mismatch vector for a j-th satellite channel, 

8 Wj represents a weight coefficient for each j-th satellite chaimel, and N represents the number of 

9 satellite channels. 

1 5. The method of claim 1 wherein step c) of generating an averaged estimate vector A* 

2 further comprises the step of: 

3 calculating said average estimate vector A4^ as A^= H ' G ' A$ c where H is a matrix of 

4 directional cosines, A* c represents said phase mismatch vector, and G is a matrix defined by 

5 G=(H^R'^H)'^H^R'^ where R is a covariance matrix of phase mismatches. 

1 6. The method of claim 1 wherein said step e) of generating a vector of controlling 

2 signals further comprises the step of: 

3 generating said vector of controlling signals in combined loop filters of each of said 

4 plurality of channels. 

1 7. The method of claim 6 where each of said combined loop filters generates a control 

2 signal ( Ucj) for a j-th channel according to the following: 

3 i 

4 Uij(i)=ETSj(k) 
5 
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6 U2j(i)=/J6j(i) + Uij(i) 

7 i 

8 U3j (i)= EU2j(k) 

10 Ucj (i)= A^j(i) + a6j (i) + U3j (i) + Vp^(i) 

11 

12 where Vprj(i) is a predicted value computed on the basis of a priori data, and a, jS, 7 are 

1 3 coefficients of the loop filter. 

1 8. The method of claim 7 wherein said coefficients are time dependent and form a 

2 Kalman-type filter. 

1 9. The method of claim 7 wherein said coefficients are constants and define a third-order 

2 channel loop. 

1 10. The method of claim 1 wherein said step g) of detecting anomalous phase 

2 measurements fiirther comprises the steps of: 

3 h) comparing residuals of said residual vector with a threshold; and 

4 i) determining that a particular channel has an anomalous phase measurement if a 

5 residual associated with said particular channel exceeds said threshold. 

1 11. The method of claim 10 wherein said step h) of comparing fiirther comprises the step 

2 of comparing the absolute value of said residuals to said threshold, wherein said threshold is in 

3 the range of 0.2 to 0.25 cycles. 

1 12. The method of claim 10 wherein said step g) of detecting anomalous phase 

2 measurements fiirther comprises the steps of: 

3 j) eliminating said particular channel firom subsequent iterations of step a); and 

4 repeating steps b) through j) until no residual exceeds said threshold or until a threshold 

5 number of said channels remain. 
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1 13. The method of claim 1 wherein said step g) of detecting anomalous phase 

2 measurements further comprises the steps of: 

3 h) generating a weighted sum of residual squares for said channels; 

4 i) comparing said weighted sum of residual squares to a threshold; and 

5 j) determining that an anomalous phase measurement exists if said weighted sum of 

6 residual squares exceeds said threshold. 

1 14. The method of claim 13 wherein said step g) of detecting anomalous phase 

2 measurements further comprises the step of: 

3 k) determining that a particular channel has an anomalous phase measurement if a 

4 residual associated with said particular channel is the maximum residual in said residual vector. 

1 15. The method of claim 14 wherein said step g) of detecting anomalous phase 

2 measurements further comprises the steps of: 

3 1) eliminating said particular channel from subsequent iterations of step a); and 

4 repeating steps b) through 1) until said weighted sum of residual squares does not exceed 

5 said threshold or until a threshold number of said channels remain. 

1 16. The mettiod of claim 13 wherein said step g) of detecting anomalous phase 

2 measurements further comprises the steps of: 

3 k) eliminating x charmels from subsequent iterations of step a), varying the channels 

4 eliminated for each iteration of step k); and 

5 repeating steps b) through k) until either said weighted sum of residual squares does not 

6 exceed said threshold or until x reaches a threshold, incrementing x by one after all combinations 

7 of channels for a given x have been eliminated during an iteration of step k). 

1 17. The method of claim 1 wherein said steps are performed in one of said satellite 

2 receivers. 

1 1 8. The method of claim 1 7 wherein said one satellite receiver receives phase 

2 measurements from the other satellite receiver. 
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1 19. The method of claim 1 further comprising the step of: 

2 h) using information about said anomalous phase measurements detected in step g) in a 

3 navigation task. 

1 20. The method of claim 19 further comprising tiie step of: 

2 i) generating a cycle slip estimate for a channel determined to have an anomalous phase 

3 measurement and using said cycle slip estimate in said navigation task. 

1 21 . An apparatus for detecting anomalous phase measurements in a satellite differential ^ 

2 navigation system in which a first satellite receiver and a second satellite receiver each compute 

3 phase measurements on a plurality of satellite channels, said apparatus comprising: 

4 a combined phase difference generator for generating a combined phase difference 

5 vector; 

6 means for generating a phase mismatch vector representing the difference between said 

7 combined phase difference vector and an estimated combined phase difference vector; 

8 an integrated discriminator for generating an averaged estimate vector representing the 

9 averaged estimate of said phase mismatch vector over said channels; 

1 0 means for generating a residual vector representing the difference between said phase 

1 1 mismatch vector and said averaged estimate vector; 

1 2 at least one joint loop filter for generating a vector of controlling signals by linear 

1 3 transformation of said residual vector and said averaged estimate vector; 

14 an accumulator for generating said estimated combined phase difference vector by 

15 successively storing components of said vector of controlling signals for each of said channels; 

16 and 

1 7 a residuals analyzer for detecting anomalous phase measurements by analyzing said 

18 residual vector. 

1 22. The apparatus of claim 21 wherein said combined phase difference vector * c 

2 comprises phase differences f*c/..*C;..*CA^) for a j-th satellite channel, and wherein said phase 

3 difference generator generates said combined phase difference vector *cas *c = *"-*'^ 

4 where vector ♦ ® comprises the full phases measured by one of said satellite receivers for each j- 
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5 th satellite channel (^bi -^ bj- -^bn) and vector * comprises the full phases measured by the 

6 other satellite receiver for each j-th satellite channel j?/. /w )■ 

1 23. The apparatus of claim 21 wherein said first and second satellite receivers are dual 

2 frequency (fl and fZ) receivers and wherein said combined phase difference vector ♦ c 

3 comprises phase differences (^ci^.^cj^.^cn) for a j-th satellite channel, and wherein said phase 

4 difference generator generates said phase differences icj of said phase difference vector * c as 

5 $cj= (4>^j - - - ^>^J)(f2 / fl) 

6 where 4> 2j and 4> ij represent the full phase measured by one of the satellite receivers at each j- 

7 th satellite channel at frequencies f2 and fl respectively, and 

8 where and $^ij represent the full phase measured by the other of the satellite receivers at 

9 each j-th satellite channel at frequencies f2 and fl respectively. 



1 24. The apparatus of claim 21 wherein said integrated discriminator generates said 

2 averaged estimate vector as 

3 N N 

4 A* = ( 2:wj-A*Cj)/( Ewj) 
5 

6 where A*Cj represents the components of said phase mismatch vector for a j-th satellite channel, 

7 Wj represents a weight coefficient for each j-th satellite channel, and N represents the number of 

8 satellite channels. 

1 25. The apparatus of claim 21 wherein said integrated discriminator generates said 

2 averaged estimate vector A^ as A4^= H * G ' A* c where H is a matrix of directional cosines, 

3 A* c represents said phase mismatch vector, and G is a matrix defined by 

4 G=(H^R'*H)"*H^R'^ where R is a covariance matrix of phase mismatches. 

1 26. The apparatus of claim 21 wherein said joint loop filter generates said vector of 

2 controlling signals in combined loop filters of each of said plurality of channels. 

1 27. The apparatus of claim 26 wherein each of said at least one combined loop filters 

2 generates a control signal ( Ucj) for a j-th channel according to the following: 
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3 i 

4 U,j(i)=E7W 

5 ^ 

6 U2j(i)=|85j(i) + U,j(i) 

7 i 

8 U3j (0= SU2j (k) 

10 Ucj (i)= A^j(i) + c^j (i) + U3j (i) + Vp^(i) 

11 

12 where Vprj(i) is a predicted value computed on the basis of a priori data, and a, /3, 7 are 

13 coefficients of the loop filter. 

1 28. The apparatus of claim 27 wherein said coefficients are time dependent and form a 

2 Kalman-type filter. 

1 29. The apparatus of claim 27 wherein said coefficients are constants and define a third- 

2 order chaimel loop. 

1 30. The apparatus of claim 21 wherein said residuals analyzer detects anomalous phase 

2 measurements by comparing residuals of said residual vector with a threshold and determining 

3 that a particular channel has an anomalous phase measurement if a residual associated with said 

4 particular channel exceeds said threshold, 

1 31. The apparatus of claim 2 1 wherein said residuals analyzer detects anomalous phase 

2 measurements by generating a weighted sum of residual squares for said channels, comparing 

3 said weighted sum of residual squares to a threshold, and determining that an anomalous phase 

4 measurement exists if said weighted sum of residual squares exceeds said threshold. 

1 32, The apparatus of claim 3 1 wherein said residuals analyzer further detects anomalous 

2 phase measurements by determining that a particular chaimel has an anomalous phase 

3 measurement if a residual associated with said particular channel is the maximum residual in 

4 said residual vector. 
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1 33. A method for detecting anomalous phase measurements in a satellite differentia^ 

2 navigation system in which a first satellite receiver and a second satellite receiver each compute 

3 phase measurements on a plurality of satellite channels, said method comprising the steps, 

4 performed for each of a plurality of iterations, of: 

5 a) generating a combined phase difference vector; 

6 b) generating an increment vector representing the difference between said combined 

7 phase difference vector and a combined phase difference vector of a preceding iteration, 

8 c) generating an averaged increment estimate vector representing the averaged estimate 

9 of said increment vector; 

10 d) generating an incremental residual vector representing the difference between said 

1 1 increment vector and said averaged increment estimate vector; 

12 e) generating an integrated residual vector from said incremental residual vector; and 

13 f) detecting anomalous phase measurements by analyzing said integrated residual vector. 

1 34. The method of claim 33 wherein step a) of generating a combined phase difference 

2 vector * c(i) further comprises the step of: 

3 calculating said phase difference vector ♦ c(i) as * c(i)= * "(i) - * '^(i) - * "^(i) where 

4 vector * "(i) comprises the fiill phases measured by one of said satellite receivers for each i-th 

5 iteration and vector * '^(i) comprises the full phases measured by the other satellite receiver for 

6 each i-th iteration and 4> ^"^(i) is a prediction of satellite movement. 

1 35. The method of claim 33 wherein said step c) of generating an averaged increment 

2 estimate vector further comprises the step of: 

A 

3 calculating said average estimate vector as 

A 

4 A<I» (i) =H(i) G(i) A*c(i), 

A 

5 where (i) represents the components of said average estimate vector for an i-th iteration, 

6 H is a matrix of directional cosines, G is a matrix defined by G=(H^R'^H)"^H^R'^ where R is a 

7 covariance matrix of phase mismatches, and A* c(i) =♦ c(i) - 4> c(i -1) where $ c(i) is a vector 

8 of combined phase difference at each i-th iteration. 
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1 36. The method of claim 33 wherein said step e) of generating an integrated residual 

2 vector further comprises the step of generating integrated residuals as 

3 5(i) = A-6(i-l) + A6(i) 

4 for each /-th iteration, where i4 is a coefficient in the range 0.995 . . . 0.999. 

1 37. The method of claim 33 wherein said step f) of detecting anomalous phase 

2 measurements further comprises the steps of: 

3 g) generating a weighted sum of integrated residual squares for said channels; 

4 h) comparing said weighted sum of integrated residual squares to a threshold; and 

5 i) determining that an anomalous phase measurement exists if said weighted sum of 

6 integrated residual squares exceeds said threshold. 

1 38. The method of claim 37 wherein said step f) of detecting anomalous phase 

2 measurements further comprises the steps of: 

3 j) eliminating x channels from subsequent iterations of step a), varying the channels 

4 eliminated for each iteration of step j); and 

5 repeating steps b) through j) until either said weighted sum of integrated residual squares 

6 does not exceed said threshold or until jc reaches a threshold, incrementing x by one after all 

7 combinations of channels for a given jc have been eliminated during an iteration of step j). 

1 39. The method of claim 33 wherein said steps are performed in one of said satellite 

2 receivers. 

1 40. The method of claim 39 wherein said one satellite receiver receives phase 

2 measurements from the other satellite receiver. 

1 41. The method of claim 33 further comprising the step of: 

2 g) using information about said anomalous phase measurements detected in step f) in a 

3 navigation task. 

1 42. The method of claim 41 further comprising the step of: 
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2 h) generating a cycle slip estimate for a channel determined to have an anomalous phase 

3 measurement and using said cycle slip estimate in said navigation task. 

1 43. The method of claim 33 further comprising the steps of: 

2 g) receiving channel indicator alarms from a channel indicator, said channel indicator 

3 alarms marking channels with anomalous phase measurements; 

4 h) eliminating channels marked with channel indicator alarms from subsequent iterations 

5 of step a); and 

6 i) repeating steps a) through f). 

1 44. The method of claim 43 wherein said step f) of detecting anomalous phase 

2 measurements further comprises the steps of: 

3 j) generating a weighted sum of integrated residual squares for said channels; 

4 k) comparing said weighted sum of integrated residual squares to a threshold; and 

5 1) determining that an anomalous phase measurement exists if said weighted sum of 



6 integrated residual squares exceeds said threshold. 

y 

1 45. An apparatus for detecting anomalous phase measurements in a satellite differential 

2 navigation system in which a first satellite receiver and a second satellite receiver each compute 

3 phase measurements on a plurality of satellite channels, said apparatus comprising: 

4 a) a combined phase difference generator for generating a combined phase difference 

5 vector; 

6 b) means for generating an increment vector representing the difference between said 

7 combined phase difference vector and a combined phase difference vector of a preceding 



8 measurement, 

9 c) an integrated converter for generating an averaged increment estimate vector 

1 0 representing the averaged estimate of said increment vector; 

1 1 d) means for generating an incremental residual vector representing the difference 

12 between said increment vector and said averaged increment estimate vector; 

13 e) a digital filter for generating an integrated residual vector from said incremental 

14 residual vector; and 
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15 f) a residuals analyzer for detecting anomalous phase measurements by analyzing said 

1 6 integrated residual vector. 

1 46. The apparatus of claim 45 wherein said combined phase difference generator 

2 generates said combined phase difference vector * c(i) as * c(i)= * ®(i) - * ^(i) - * ^"^(1) where 

3 vector * ^(i) comprises the full phases measured by one of said satellite receivers for each i-th 

4 iteration and vector * '^(i) comprises the full phases measured by the other satellite receiver for 

5 each i-th iteration and * ^'^(i) is a prediction of satellite movement. 

1 47. The apparatus of claim 45 wherein said integrated converter generates said averaged 

A 

2 increment estimate vector as 

A 

3 A<I> (i) =H(i) G(i) A*c(i), 

A 

4 where (i) represents the components of said average estimate vector for an i-th iteration, 

5 H is a matrix of directional cosines, G is a matrix defined by G=(H^R"*H)''H^R"^ where R is a 

6 covariance matrix of phase mismatches, and A* c(i) c(i) - * c(i -1) where * c(i) is a vector 

7 of combined phase difference at each i-th iteration. 

1 48. The apparatus of claim 45 wherein said residuals analyzer detects anomalous phase 

2 measurements by generating a weighted sum of integrated residual squares for said channels, 

3 comparing said weighted sum of integrated residual squares to a threshold, and determining that 

4 an anomalous phase measurement exists if said weighted sum of integrated residual squares 

5 exceeds said threshold. 



1 49. The apparatus of claim 45 wherein said digital filter generates integrated residuals 

2 £(i) as 

3 6(i) = A-5(i-l) + A6(i) 

4 for each of an /-th iteration, where /I is a coefficient in the range 0.995 . . . 0.999. 

1 50. A method for detecting anomalous phase measurements in a satellite differential 

2 navigation system in which a first satellite receiver and a second satellite receiver each compute 
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3 phase measurements on a plurality of satellite channels, said method comprising the steps, 

4 performed for each of a plurality of iterations, of: 

5 a) generating a combined phase difference vector; 

6 b) generating an increment vector representing the difference between said combined 

7 phase difference vector and a combined phase difference vector of an initial measurement at an 

8 initial time, 

9 c) generating a corrected increment vector using a cycle slip correction estimate 

1 0 generated in a preceding iteration; 

11 d) analyzing channel indicator alarms; 

12 e) generating an averaged increment estimate vector representing the averaged estimate 

13 of said corrected increment vector using channels not associated with a channel indicator alarm; 

14 f) generating a residual vector representing the difference between the corrected 

1 5 increment vector and the averaged increment estimate vector; and 

16 g) generating said cycle slip correction estimate using said residual vector. 

1 51. The method of claim 50 wherein said step a) of generating a combined phase 

2 difference vector further comprises the step of: 

3 including a prediction of phase difference in said difference. 

1 52. The method of claim 50 wherein said initial time is reset periodically. 

1 53. The method of claim 50 further comprising the step of: 

2 h) using said cycle slip correction estimate in a navigation task. 

1 54. An apparatus for detecting anomalous phase measurements in a satellite differential 

2 navigation system in which a first satellite receiver and a second satellite receiver each compute 

3 phase measurements on a plurality of satellite charmels, said apparatus comprising: 

4 a) a combined phase difference generator for generating a combined phase difference 

5 vector; 

6 b) means for generating an increment vector representing the difference between said 

7 combined phase difference vector and a combined phase difference vector of an initial 

8 measurement at an initial time; 
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9 c) a cycle slip correction unit for generating a cycle slip correction estimate; 

1 0 d) means for generating a corrected increment vector using a cycle slip correction 

1 1 estimate generated in a preceding measurement; 

12 d) a channel indicator analyzer for analyzing channel indicator alarms; 

13 e) an integrated converter for generating an averaged increment estimate vector 

14 representing the averaged estimate of said corrected increment vector using channels not 

1 5 associated with a channel indicator alarm; 

16 f) means for generating a residual vector representing the difference between the 

17 corrected increment vector and the averaged increment estimate vector; and 

18 g) a correction unit for generating said cycle slip correction estimates using said residual 

19 vector. 

1 55. The apparatus of claim 54 wherein said combined phase difference generator further 

2 generates said combined phase difference vector by including a prediction of phase difference in 

3 said difference. 

1 56. The apparatus of claim 54 wherein said initial time is reset periodically. 



Page 38 of 39 



